Percolative conductivity and critical exponents in mixed-valent manganites 
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Recent experiments have shown that some colossal magnetoresistance (CMR) materials exhibit a 
percolation transition. The conductivity exponent varies substantially with or without an external 
magnetic field. This finding prompted us to carry out theoretical studies of percolation transition 
in CMR systems. We find that the percolation transition coincides with the magnetic transition 
and this causes a large effect of a magnetic field on the percolation transition. Using real-space- 
renormalization method and numerical calculations for two-dimensional (2D) and three-dimensional 
(3D) models, we obtain the conductivity exponent t to be 5.3 (3D) and 3.3 (2D) without a magnetic 
field, and 1.7 (3D) and 1.4 (2D) with a magnetic field. 

PACS numbers; 75.90.-fw, 75.45+j, and 71.27.+a 



Colossal magnetoresistance (CMR), an unusual large 
change of resistivity in the presence of a magnetic field, 
has been extensively studied in ferromagnetic perovskite 
manganites. It is well known that the electronic 
phase diagrams of CMR materials are very complex. 
There are various ordering states and phase transitions 
as the carrier concentration is varied. For example, 
Lai_a;Caa;Mn03, a typical double exchange ferromag- 
netic metal when x=0.33, becomes a charge ordered in- 
sulator with a specific type of electronic orbital and 
magnetic orderings when x=l/2. |^ Recently, it was 
demonstrated that Lai_a;_j,Prj,Caa;Mn03 (x=3/8) sys- 
tem, where Pr is chosen to vary the chemical pres- 
sure, may be electronically phase separated into a sub- 
micrometer-scale mixture of insulating regions and fer- 
romagnetic (FM) metallic regions. Experimental 
findings of FM clusters and phase separation in CMR 
materials have also been reported in early studies. 
1^ Electron diffraction and dark-field imaging on the 
Lai_a;_j,PryCa3;Mn03 samples indicate that the insulat- 
ing region is a x=l/2 charge ordered phase. This is not a 
charge congregation type of phase separation, which was 
observed in slightly doped antiferromagnetic manganites 
1^ and was extensively discussed. [Q CMR effect was ob- 
served in different samples with 0.275 < y < 0.41, and 
was explained by percolative transport through the ferro- 
magnetic domains. According to the percolation theory, 
the conductivity a oc (p — po)*, where p is the concen- 
tration of the metallic phase and po is its critical value. 
In Ref . 1^ , the exponent was studied in the presence and 
absence of an external magnetic field and two values are 
substantially different. Thus, the experiment shows that 
the percolative transport in the CMR systems depends 
sensitively on the relative spin orientation of adjacent 
ferromagnetic domains which is controlled by an applied 
magnetic field. The goal of this work is to develop a per- 
colation theory which takes into account the magnetic 
transition in the CMR materials. In particular, we would 
like to understand what causes the exponent to be so dif- 



ferent with or without an external magnetic field. 

Based on phenomenological considerations, we study 
two- or three-dimensional (2D or 3D) lattice percolation 
models. The conductivity and its critical exponent are 
calculated by means of real space renormalization and 
numerical methods. The system consists of three types 
of lattice sites. Each site has spin Si{= ±1,0) where 
i denotes the index of a site. This model is similar to 
the site-diluted spin system used in Ref. Si = 

implies that the site is empty, meaning occupied by the 
x=l/2 charge ordered phase. Si = ±1 means the site 
is occupied by the ferromagnetic metallic phase with up 
and down magnetizations. The Hamiltonian of the 
spin interaction is Ising-like and is written as 



-jj2s. 



s, + J2h-s. 



(1) 



Here denotes a pair of the nearest neighbor sites, J 
is the interaction energy, and H is the strength of the 
external magnetic field. We consider a ferromagnetic 
interaction, i.e., J > 0. In order to investigate trans- 
port properties, we assume that the local conductivity 
between two nearest occupied sites is either for an- 
tiparallel spins (ti and J,t) or 1 for parallel spins (IT and 
ID and is zero if one or both sites are empty. Thus, the 
conductivity between neighboring sites can be expressed 
as follows 



1, Si ■ Sj — 1; 
"10, otherwise. 



(2) 



First, we use the standard Monte-Carlo method to 
study the magnetic properties of the model. Knowing 
the spin structure is necessary for the transport studies 
since the local conductivity depends on the spins through 
Eq.(2). The conductance calculation is performed on a 
set of spin configurations produced by the Markov chain. 
In the Markov chain every spin configuration is gen- 
erated from the previous one by using the probability 
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^-l3AHs ^^g-fiAHs _l_ j^-j^ where AHg is the energy differ- 
ence between these two configurations and P = l/kgT. 
pTf The calculations are carried out on finite square and 
cubic lattices for 2D and 3D systems and the periodic 
boundary condition is adopted to eliminate the boundary 
effects. The conductance G for every spin configuration 
in the Markov chain is obtained by calculating the total 
conductance of the resistor network. In the resis- 
tor network, the local conductivity between neighboring 
sites is determined by Eq.(2). The total conductance is 
calculated for many samples to obtain the average con- 
ductance. 

When the charge order (CO) phase is dominant in 
some of the CMR materials, most sites are empty ac- 
cording to our model. Thus, most of the conductivities 
between neighboring sites are zero, corresponding to the 
low concentration limit {p 0) in the percolation model, 
where p is the probability of non-zero local conductiv- 
ity. The small FM islands is well separated by the CO 
phase. Because the spin correlations between FM blocks 
are cut off by the CO phase, the spin orientation for each 
FM block is random, either up or down. Therefore, the 
spontaneous magnetization m will be zero at any tem- 
perature T. As p increases, m continues to be zero until 
p reaches pa at which the first infinite FM cluster ap- 
pears. If p > poi a- finite spontaneous magnetization ap- 
pears for T < Tc{p) with Tcij)) the critical temperature 
at concentration p. In Figs. 1(a) and 1(b) we plot the 
normalized magnetization m as a function of tempera- 
ture T with different concentration probabilities p for 2D 
and 3D systems, respectively, calculated using the Monte 
Carlo method. The interaction strength J is set to unity. 
It can be seen that when the temperature increases m is 
reduced and reaches zero at Tdp). From magnetization 
data with different size samples and through finite-size 
scaling, we can determine the critical temperature. In 
Figs. 2(a) and 2(b) we plot the relation between the criti- 
cal temperature Tc{p) and the concentration p. As shown 
in these plots, Tdp) will approach the Curie temperature 
Tc for the regular Ising model as p approaches to 1. We 
find that Tc determined from Fig. 2 is 2.3 J or 4.4 J for 2D 
or 3D system. Both of them are in good agreement with 
the values of the regular Ising model. |l^ 

From the spin configurations of the system, we can cal- 
culate the conductance by using the local conductivity 
defined in Eq. (^). In a paramagnetic phase, there is no 
infinite cluster with spins pointing in the same direction. 
Hence, according to Eq. (y) there is no conducting path 
throughout the sample. On the other hand, the first con- 
ducting path appears simultaneously when the magneti- 
zation starts becoming non-zero. Thus, the phase tran- 
sition from the FM phase to the paramagnetic phase is 
accompanied by the metal-insulator (MI) transition in the 
conductance. This implies that at zero temperature the 
MI transition occurs at the critical concentration of the 
percolation threshold pq. Near the critical point the av- 
eraged conductance G can be expressed as G ~ {p — poY, 
where t is the the conductance critical exponent for the 
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FIG. 1. (a)The normalized spontaneous magnetization m 
as a function of temperature T for different concentration p 
in the 2D model. The unit of T is set by the interaction 
strength J/ks- From top to bottom, the concentrations are 
p — 1.0, 0.9, 0.8, 0.7, 0.6 and 0.5 The calculation is performed 
on 1000 samples whose size is 100 x 100. (b) m as a function 
of T in the 3D model. The calculation is performed on 200 
samples whose size is 10 x 10 x 10. From top to bottom, the 
concentrations are p = 1.0, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3 and 
0.2. 
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FIG. 2. (a) The critical temperature as a function of 
the concentration p in 2D. (b) The critical temperature Tc as 
a function of the concentration p in 3D. 




FIG. 3. To illustration of renormalization in 2d triangle 
lattice 



transition. 

We now discuss the real space renormalization method. 
We start with a 2D system and later extend to a more 
complicated 3D system. The essential physics does not 
depend on dimensionality, although numerical numbers 
do. Consider a triangular lattice in a 2D plane. The 
choice of a triangular lattice is for the convenience of 
rescaling of the system in the real space renormaliza- 
tion procedure. By enlarging the system by factor VS 
and grouping every three sites into a " supersite" , the 
number of the supersites is the same as the number of 
sites in the original system (see Fig. 3). In Fig. 3, the 
thin lines are for the original system and the thick lines 
are for the rescaled system. A supersite is regarded as 
empty if the majority of the three sites are empty. Thus 
the concentration of the new system can be expressed 
as p' = p^ + 3p'^{l - p)(Ref. i|l|l). Near the critical 
point Po — 0.5, on the metallic side, G can be written as 
G = Go(p— po)*j where the constant Gq is proportional to 
the conductance of the unit cell. For the enlarged system 
this relation becomes G" = Gq{p' —poY- In a 2D system 
the conductance is independent of the system size and we 
have G = G'. Thus, t = ln(ga)/in(^££a). if the spin 
degree of freedom is frozen, the conductance of elemen- 
tary cell is inversely proportional to the size of the cell, 
hence, t = ln(\/3)/ ln(|) = 1.35. This value is very close 
to the exact value of i = 4/3 in the standard 2D per- 
colation model. ||l^ When taking into account the spin 
degrees of freedom, as we discussed above, the pq point is 
also the critical point for the spontaneous magnetization. 
The conductance Go (or Gq) is now associated with spin 



configurations. The spin dependence comes from the fact 
that Go is proportional to the average of conductivity be- 
tween two nearest neighbor sites which is spin dependent, 



Go 



cos(- 
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^—f3mS cos 9i ^—{5mS cos 9j 



cos(^) 



(3) 



The conductivity expression between nearest neigh- 
bor spin try ~ cos(^) comes from the double-exchange 
model where Oij is the angular difference between spins 
Si and Sj , which satisfies |15 



cos 9ij — cos 6i cos 6j + sin 9i sin 9j cos (0i — ) (4) 

{9i,(j)i) denotes the orientation of spin Si. A is a nor- 
malization constant. From the above equation, it is easy 



to show that Go 



Because m can be written as 



m (p — po), |13[ we finally get the conductance critical 
exponent in 2D triangular lattice as ~ 1.35 -f 2 = 3.35. 

In the 3D case, we consider the normal cubic lattice. 
The elementary vectors of the enlarged lattice is just 
two times of those of the original one {cx' ,ey' ,ez'} = 

{2e-;, 2ey, 2e J. In this case p' = En=4 C8>"(1 - p)"*"", 
from which we conclude that po = 0.395 and t = 1.7 
without spin effects. After the spin degrees of freedom 
are considered, the formula Go ~ rri^ is still satisfied. But 
in 3D case m has m ^ (p — Po)"'^ ™- So the critical 
exponent t' ~ 1.7 + 2 x 1.79 ~ 5.3. 

In Fig. 4 we show the numerical results of 2D con- 
ductance G. 3D calculation has not been done because 
of computational limit. The points with steepest drop 
are defined as the critical points for the metal-insulator 
transition. These points are consistent with the magnetic 
critical points (see Fig.l), as we discussed before. In the 
inset of Fig. 4 we show the ln(G) versus \n{p — po). They 
exhibit a linear dependence and the slope of the curve, 
corresponding to the exponent i, is roughly 3.5. This 
value is in a good agreement with the estimation from 
the renormalization group consideration although it is 
obtained from a different type of lattice structure. In 
Table I we list the results of this work and compare them 
with the previous studies of the standard percolation the- 
ory jlZj and experimental exponents in CMR materials. 

The experimental exponent (or t'^) was obtained 
in the absence (or presence) of an external magnetic field. 

Before summary, we would like to make couple com- 
ments, (i) In order to calculate the conductivity ex- 
ponent, we have to determine the dependence of Go 
(see Eq.(3)) on the magnetization to, and the double- 
exchange model is used to achieve that goal. However, 
the double-exchange may not be the origin of the ferro- 
magnctism in doped manganites, as shown in a recent 
work, mm This might explain some of the discrepancy 
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FIG. 4. The normalized conductivity G as a function of 
temperature T in 2D model. From top to bottom, the con- 
centrations are p — 1.0,0.9,0.7,0.6 and 0.5. Inset: logG vs. 
log(p — po) shows a linear behavior with a slope ~ 3.5 

TABLE I. The critical exponent t for 2d and 3d cases. Pc2 
and Pc3 are percolation threshold for 2d and 3d, respectively. 
td and t'd are the critical exponents for d-dimension model 
without or with spin effects. Previous results are from Ref.l4. 
Experiment I is from Ref.3 and Experiment II is from Ref.4. 



Previous results 
Experiment I 
Experiment II 
Renormalization method 
Numerical results 



Pc2 

0.5 



0.5 
0.51 



Pc3 

0.31 



0.37 
0.32 



t2 

1.33 



1.35 
1.30 



t3 

1.9 
2.6 

1.7 



t'. 



3.3 
3.5 



t's 

6.1 
2.6 
5.3 



between the experimental and our theoretical values in 
the conductivity exponent. We should also mention that 
the values from the two experimental groups are not in 
good agreement with each other as seen in Table I. (ii) 
We have neglected the quantum effects in this work. This 
might be justified because Tc is relatively high in these 
samples. Develop a semi-classical transport theory in 
this problem is a difficult task because of the finite phase 
coherence length. There are attempts |17 of using the 
semi-classical theory to understand the two-dimensional 
metal-insulator transition. 

In summary, we argue that the percolation threshold 
corresponds not only to the appearance of an infinite 
metal cluster, but also to the phase transition from PM to 
FM. This coincidence of two phase transitions renormal- 
izes the critical exponents. The conductivity exponent 
has been obtained using real space renormalization and 
numerical calculations. The exponent is found to be quite 
different whether the magnetic transition is considered. 
This finding explains the large exponent discrepancy in 
some of the CMR materials in the presence or absence of 



an external magnetic field. 
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